Efficient Quasiparticle Evacuation in Superconducting Devices 
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We have studied the diffusion of excess quasiparticles in a current-biased superconductor strip in 
proximity to a metallic trap junction. In particular, we have measured accurately the superconductor 
temperature at a near-gap injection voltage. By analyzing our data quantitatively, we provide a full 
description of the spatial distribution of excess quasiparticles in the superconductor. We show that 
a metallic trap junction contributes significantly to the evacuation of excess quasiparticles. 



In a normal metal - insulator - superconductor (N-I-S) 
junction, charge transport is mainly governed by quasi- 
particles The presence of the superconducting en- 
ergy gap A induces an energy selectivity of quasiparticles 
tunneling out of the normal metal 0, Q ■ The quasiparti- 
cle tunnel current is thus accompanied by a heat trans- 
fer from the normal metal to the superconductor that is 
maximum at a voltage bias just below the superconduct- 
ing gap (V < A/e). For a double junction geometry (S- 
FN-FS), electrons in the normal metal can typically cool 
from 300 mK down to about 100 mK @, 0, llj- However, 
in all experiments so far the electronic cooling is less effi- 
cient than expected 0, [j| . It has been proposed that this 
inefficiency is mostly linked to the injected quasiparticles 
accumulating near the tunnel junction area. This out-of- 
equilibrium electronic population, injected at an energy 
above the superconductor energy gap A, relaxes by slow 
processes such as recombination and pair-breaking pro- 
cesses. The accumulation of quasiparticles is aggravated 
in sub-micron devices, where the relaxation processes are 
restricted by the physical dimensions of the device, lead- 
ing to an enhanced density of quasiparticles close to the 
injection point. These quasiparticles can thereafter tun- 
nel back into the normal metal Q , genera ting a parasitic 
power proportional to the bias current '7, 12]. The same 
phenomenon is relevant to other superconductor-based 
devices such as qubits [gj], single electron transistors [i| 
and low temperature detectors [lOj . 

In hybrid superconducting devices fabricated by mul- 
tiple angle evaporation, a normal metal strip in tunnel 
contact with the superconducting electrode acts as a trap 
for excess near-gap quasiparticles, which removes them 
from the superconductor. This mechanism is usually not 
fully efficient due to the tunnel barrier between the nor- 
mal metal and the superconductor FLj. A detailed the- 
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FIG. 1: (a) Schematic of the sample design with a trap junc- 
tion on each S-strip of the superconducting microcooler. The 
curve shows the spatial profile of the excess quasiparticles 
SN qp along the S-strip. The two injectors are located at ai = 
5 fim or a.2 — 21 fim; the detector (cooler junction island) is 
at x = fim and the length of the S-strip is Ls = 27 /im. (b) 
Calculated spatial decay of excess quasiparticles density 5N qp 
along the S-strip. The two injectors are biased at eV/A=4. 



of the diffusion of out-of-equilibrium quasiparticles in a 
superconducting strip covered with a trap junction. A 
N-metal is used to inject quasiparticles in the S-strip. 
The local superconductor temperature is inferred from 
the heating of the central N-island of a S-FN-I-S junction. 
We quantitatively compare our experimental data with a 
recently discussed theoretical model 0. 

We have used a S-I-N-I-S cooler device with a trap 
geometry similar to the one studied in Ref. [13j, see 
Fig. 1. These devices are fabricated using electron beam 



ory of non-equilibrium phenomena in a superconductor 
in contact with normal metal traps has been developed 
fi"3 |. However, a quantitative comparison between exper- 
iments and theoretical predictions is so far still missing. 

In this Fetter, we present an experimental investigation 



lithography, two-angle shadow evaporation and lift-off on 
a silicon substrate having 500 nm thick Si02 on it. The 
central normal metal Cu electrode is 0.3 /im wide, 0.05 
/im thick and 4 /im long. The 27 /im long symmetric S- 
strips of Al are then partially covered, through a tunnel 
barrier, by a Cu strip acting as a trap junction. At their 
extremity, the S-trips are connected to a contact pad act- 
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ing as a reservoir. In addition to the cooler island, we 
added two normal metal Cu tunnel injector junctions of 
area around 0.09 /mi 2 on one S-strip. Injector 1 and 2 are 
at a distance of a = 5 (im and 21 /im respectively from 
the central Cu island. The Al tunnel barrier is assumed 
to be identical in the cooler, probe and trap junctions, 
since they have similar specific tunnel resistance. The 
normal state resistance of the double-junction S-I-N-I-S 
cooler is 1.9 kf2. The normal state resistance of N-I-S 
injector junctions 1 and 2 are respectively 2.5 kS7 and 2.3 
kil. The diffusion coefficient of the Al S-strip film was 
measured at 4.2 K to be 30 cm 2 /s. 

N-I-S tunnel junctions are known to enable controlled 
quasiparticle injection in a superconductor 14, TBI [2l|. 



The tunnel current through a N-I-S junction is given by: 



I(V) 



n s (E)[f N {E ~ eV) - f N (E + eV)]dE 

(1) 

where Rn is the normal state resistance, /n is the elec- 
tron energy distribution in the normal metal and rig is 
the normalized density of states in the superconductor. 
In a superconducting wire undergoing quasiparticle in- 
jection, the superconductor gap A(Tg) is suppressed lo- 
cally. As this gap can be extracted from a N-I-S junc- 
tion current-voltage characteristic, such a junction can 
be used for quasiparticles detection. Usually, an effec- 
tive superconductor temperature Tg is inferred from the 
superconductor gap-temperature dependence. Fig. 2(a) 
displays the differential conductance dl/dV of a N-I-S 
probe junction (similar to an injector in Fig. 1(a)) lo- 
cated on a S-strip at different cooler bias voltages. At 
high injection, the gap suppression appears clearly, and 
enables a good determination of the superconductor ef- 
fective temperature. This a ppr oach was used in numer- 
ous previous studies 14, IH 21 1. At lower injection with 
a voltage closer to the gap voltage, the tunnel charac- 
teristic becomes little sensitive to quasiparticle injection. 
For instance, in Fig. 2a the probe junction characteristic 
at 1 mV injection voltage almost overlaps the equilibrium 
characteristic (at mV). This limitation comes naturally 
from the saturation of the superconducting gap at low 
temperature Tg^; Tc, where Tq is the superconductor 
critical temperature. So far, the insensitivity of the N-I- 
S junction characteristic at low injection bias has been a 
major roadblock in investigating the decay of quasipar- 
ticles injected at energies just above the gap la. \21\. 

Instead of measuring directly the superconductor, a 
better detection sensitivity can be achieved by measuring 
the temperature of a small N-island connected to super- 
conductor through a tunnel barrier [2(| • In the absence of 
excess quasiparticles in the superconductor, the N-island 
is in thermal equilibrium with it. When the superconduc- 
tor is under injection, some of the excess quasiparticles 
population will escape by tunneling (even at zero bias) 
from the superconductor to the central N-island. The in- 



jected quasiparticles population will then reach a quasi- 
equilibrium in the N-island with a electronic temperature 
Tjyj different from the cryostat temperature T^ a ^. 
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FIG. 2: (a) Probe junction differential conductance under 
different injection bias voltage from the S-I-N-I-S cooler junc- 
tion, (b) Low bias cooler junction data (solid black lines) 
at different injector 2 voltage level compared to calculated 
isotherms (red dashed lines) as obtained from Eq. 1 with A 
= 0.22 meV at different T N . 

In this study, N-I-S junctions located on one S-strip 
of a S-I-N-I-S junction are used as quasiparticle injectors 
by current-biasing them. This leads to a spatial distribu- 
tion of the excess quasiparticles density along the S-strip 
(see Fig 1 (b)). The central N- metal in the S-I-N-I-S 
cooler geometry is used as a detector for the quasipar- 
ticle density (at x — 0) in the superconductor. In the 
N-island, the phase coherence time of about 200 ps (mea- 
sured from a weak localization experiment in a wire from 
the same material) is much shorter than the mean escape 
time from the island estimated to about 100 ns. The N- 
island electronic population is then at quasi-equilibrium. 
Its temperature Tj^ can be extracted from the zero bias 
conductance level of the S-I-N-I-S junction Further, 
the superconductor temperature is inferred from the N- 
metal temperature by considering its heat balance. As 
demonstrated below, this scheme is highly sensitive down 
to about 200 mK, where the S-I-N-I-S junction I-V be- 
comes dominated by the Andreev current [l3j]. 

Fig. 2(b) shows the differential conductance of the 
cooler junction (full black lines) at different injector-2 
bias voltages along with isotherms (dotted red lines) cal- 
culated from Eq. 1. Fig. 3(a) displays the central N- 
metal temperature extracted from the zero-bias conduc- 
tance as a function of injector bias voltage. As the injec- 
tor bias increases above the bath temperature 2b a th' ^ ne 
temperature increases, indicating that more quasi- 
particles tunnel from the S-trip to the N-island. 

In order to obtain the superconductor temperature Tg 
at the cooler edge (x = 0), we need to consider the heat 
balance in the normal metal. The heat flow across a N- 
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I-S junction with different quasiparticle distribution on 
either side of the tunnel barrier is given by: 



Pheat(T N ,Ts) 



e 2 R N 



' ' En s (E)[f N (E) - f s (E)]dE 

where fs is the energy distribution function in the su- 
perconductor at temperature Tg. It is compensated by 
electron-phonon coupling power P e - P h so that Pheat + 
P e - P h — 0- Here, we have used the usual expression 
for the electron-phonon coupling P e - P h = T*U(T^ — 
T% h ), where S = 2 nW ■ fim,- 3 ■ K~ 5 in Cu is a 
material-dependent constant and U is the metal vol- 
ume. For the normal metal phonons, the electron- 
phonon coupling power is compensated by the Kapitza 
power P K (T ph ; T bath ) = KA(T A ath - T p \), where K is an 
interface-dependent parameter and A the contact area. 
The inset of Fig. 3(a) displays the correspondence be- 
tween the superconductor temperature Tg and the cen- 
tral N- metal temperature T^. We took the fitted Kapitza 
coupling parameter value K.A = 144 pW ■ K~ A found in 
Ref. [13J in a very similar sample. The grey area shows 
the calculated temperature Tg at different Kaptiza cou- 
pling coefficient K ranging from 120 W ■ mT 2 ■ K~ A to 
infinity. The uncertainty in Tg due to uncertainty in K 
is negligible below 400 mK and at Tg at 600 mK it is 
only 15 mK. 
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FIG. 3: (a) Central N-metal electronic temperature Tjv de- 
pendence on the injector 2 bias voltage at the bath tempera- 
ture of 304 mK. Error in Tjv is smaller than the symbol size. 
Inset: calibration of the extracted superconductor tempera- 
ture Ts to the measured normal metal temperature Tat. The 
grey area shows the uncertainty in Ts for different values of 
the Kapitza coupling coefficient, (b) Corresponding extracted 
S-metal temperature Ts at x = 0. Dotted lines are fits using 



35 cm 2 /s. 



Fig. 3(b) and Fig. 4 show the extracted supercon- 
ductor temperature Ts at the cooler edge (a; = 0) for 
different injection bias, in the two injectors at three dif- 
ferent bath temperatures T^ a ^ n = 100 mK, 304 mK and 
500 mK respectively. We have succeeded in obtaining 



accurately the superconductor temperature for injection 
bias voltage close to the gap voltage. For a given injec- 
tion bias in the two injectors, the temperature Tg of the 
detector at x = is higher for a closer injector. This con- 
firms qualitatively the diffusion-based relaxation of hot 
quasiparticles in the superconductor. 

In an out-of-equilibrium superconductor, the quasi- 
particle density N qp and the phonon density of 2A en- 
ergy A^2A are coupled to each other by the well-known 
Rothwarf- Taylor (R-T) equations [lg. In a recent work 
[7(, some of us have extended the R-T model to include 
the influence of the trap junction on the quasiparticle dif- 
fusion. We considered a superconducting strip covered 
by a second normal metal separated by a tunnel barrier, 
which is in practice equivalent to a device fabricated with 
a shadow evaporation technique Q . The normalized ex- 
cess spatial quasiparticle density z(x) in the S-strip is 
given by the solution of the differential equation: 



D„ 



d 2 ^ 
' dx 2 



z z + z 2 /2 
— + — . 

T T eff 



(3) 



From this equation, one can find the quasiparticle de- 
cay length A = y/D qp T e ff /a, where r e // is the mate- 
rial dependent effective recombination time and a = 1 + 
T e ff /to (>1) is the enhancement ratio of the quasiparti- 
cle decay rate due to the presence of the trap junction. 
When the trapping effect is dominant, the quasiparticle 
decay length reduces to U D qp To. The trap characteristic 
time To describes the rate of quasiparticles escaping to 
the N-metal trap. It is defined as t = e 2 RNNN(Ep)ds, 
where Rnn is the specific resistance of the trap junction 
and ds is the thickness of S. At equilibrium, the den- 
sity of quasiparticles N qp o in the superconductor at tem- 
perature Tg (< Tq) decays exponentially: N qp o(Ts) — 
N(E F )A(nkT s /2A) 1 / 2 exp[-A/fcT s ], where N(E F ) is 
the density of states at the Fermi level. 

In our experiment, the injectors are biased just above 
the gap voltage. Here, we assume that the injected quasi- 
particles relax fast (in comparison to other recombination 
processes) to the superconductor gap energy level and 
thus can be afterwards adequately described by the cou- 
pled R-T equations. To compare our experimental result 
with the theoretical model, we have extended the de- 
scription shown in Ref. Q ■ We solved Eq. 3 numerically 
with boundary conditions so as to include the injection: 
$|| + - |f |_ = %i at x = a; the detection: || = at 
x = 0; and the finite length of the S-strip: z = at x — 
Lis into the model. The last boundary condition x = L5 
provides an additional path for the excess quasiparticles 
to thermalize in addition to the N-trap. Further, the lo- 
cally enhanced quasiparticle density N qp {x) is described 
by an equilibrium quasiparticle density N qp Q such that 
N qp (x = 0) = N qp o(T = Ts) to obtain the theoretical 
superconductor temperature Tg (l7| . 

In the fit procedure, we used the calculated values of 
T e f f and tq . The calculated effective recombination time 
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FIG. 4: Extracted superconductor temperature Ts as a func- 
tion of injector bias voltage at a cryostat temperature of 100 
mK (square dots) and 500 mK (circle dots), (a) and (b) 
correspond to an injection from injector 1 and 2 respectively. 
Dotted lines show the fits with the parameter D qp — 35 cm 2 /s. 
Red and blue dashed lines are the calculated curve with the 
same D 9P = 35 cm 2 /s and for tq = 12 ns and 00 respectively. 



for Al is T e ff = 14 /Lts 
imental value in Ref. 



181 ] which is close to the exper- 
19(. For our sample parameters, 
the calculated trap characteristic time To is equal to 0.3 
/is. The model has then only one free parameter, which 
is the quasiparticle diffusion coefficient D qp . We obtain 
a quantitative agreement between the theoretical predic- 
tions for D qp = 35 cm 2 /s and the experiment (dotted 
lines in Fig. 3(b) and Fig. 4) on the two injectors, for 
every injection voltage, and for a bath temperature of 
100 mK to 500 mK. The fit-derived value of T) qp is com- 
parable to the measured diffusion coefficient for Al at 4.2 
K, and corresponds to A = 30 /im. 

This excellent agreement demonstrates that we quanti- 
tatively understand the diffusion of quasiparticles in the 
superconductor in proximity with the metallic trap junc- 
tion. The quasiparticles relax mostly through two chan- 
nels: N-metal trap (decay length A) and absorption in 
the reservoir (decay length Lg). In our device, the decay 
length of both channels is around 30 /im, thus they act 
with a similar efficiency. Fig. 4 red and blue dashed lines 
shows the calculated superconductor temperature Tg at 
-^bath = m ^ f° r a more transparent trap junction 
To = 12 ns and 00, in parallel to the experiment and its 
fit. For To = 00, the excess quasiparticles relax at the 
reservoir x — Lg. For tq = 12 ns, the decay length is A 
= 6.5 /mi, so that trapping then dominates the quasipar- 
ticles absorption. At such transparency the influence of 
proximity effect cannot be ignored 

In conclusion, we have studied experimentally the dif- 
fusion of quasiparticles injected in a superconductor with 
an energy close to the gap voltage A in the presence 
of normal metal trap traps. Our study demonstrates 
that in such devices quasiparticle trapping competes with 
relaxation in the reservoir. This new knowledge is of 



great importance in improving the geometry of the fu- 
ture cooling devices and other superconductor, based 
low-temperature devices. 
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